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While supernova remnants (SNRs) are widely thought to be powerful cosmic-ray accelerators, 
indirect evidence comes from a small number of well-studied cases. Here we systematically de¬ 
termine the gamma-ray emission detected by the Fermi Large Area Telescope (LAT) from all 
known Galactic SNRs, disentangling them from the sea of cosmic-ray generated photons in the 
Galactic plane. Using LAT data we have characterized the 1-100 GeV emission in 279 regions 
containing SNRs, accounting for systematic uncertainties caused by source confusion and in¬ 
strumental response. We have also developed a method to explore some systematic effects on 
SNR properties caused by the modeling of the interstellar emission (lEM). The lEM contributes 
substantially to gamma-ray emission in the regions where SNRs are located. To explore the sys- 
tematics we consider different model construction methods, different model input parameters, and 
independently fit the model components to the gamma-ray data. We will describe this analysis 
method in detail. In the Eirst Fermi LAT SNR Catalog there are 30 sources classified as SNRs, 
using spatial overlap with the radio position. Eor all the remaining regions we evaluated upper 
limits on SNRs’ emission. In this work we will present a study of the aggregate characteristics of 
SNRs, such as comparisons between GeV and radio sizes as well as fluxes and spectral indexes 
and with TeV. 
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1. Introduction 

In this work we will shortly describe the procedure used in the systematic search and analysis 
of supernova remnants (SNR) in the F<?rm/-Large Area Telescope (LAT) data. The complete pro¬ 
cedure will be described in a future publication [1]. Several results of this work will be described 
and highlighted also in [2]. 

2. Analysis procedure 

Our analysis starts by considering the 274 radio SNRs described in [3] and five new SNRs dis¬ 
covered later. Using their radio position and size we analyze three years of Pass7v6 “Source”-class 
Fermi LAT data in the energy range between 1 and 100 GeV. For each of these radio SNRs we test 
their possible gamma-ray counterpart extension, location, spectral curvature and significance. The 
method used for the spectral and spatial fit is described in [4]. Since we don’t know exactly the 
background sources a priori for all these sources we developed a method that iteratively adds back¬ 
ground sources on top of the standard interstellar emission model (IBM) gal_2yearp7v6_v0.fits, 
and isotropic model iso_p7v6source.txt. As background sources we used also the pulsar described 
in the second pulsar catalog [5] and the associated sources in the LAT second source catalog [6]. 
For all this candidate SNRs we evaluated their spectral characteristics and location, if we did not 
detect any source we reported a 99% Bayesian upper limit (UL) [7] . 

If GeV and radio sizes are similar, as has been observed on an individual basis for several ex¬ 
tended SNRs (e.g. [4]), the LAT has sufficient spatial resolution to detect many SNRs as extended. 
Figure 1 shows the distribution of radio diameters from Green’s catalog. Vertical dashed lines show 
the minimum detectable extension for source with flux and index typical of those observed in this 
catalog, based on simulations using the P7V6 instrument response functions (IRFs) [4]. 

3. Detection method and source classification 

For each SNR, we character size, morphology and spectrum of any y-ray emission that may 
be coincident with the radio position reported in Green’s catalog. This was achieved by testing 
multiple hypotheses for the spatial distribution of y-ray emission: a point source and two different 
extended disk. One was obtained using only a disk source coincident in size with the radio size, 
while the second was obtained allowing the nearby background sources to be included in the final 
size of the SNR. The best fit was selected based on the global likelihoods of the fitted hypotheses 
and their numbers of degrees of freedom. To define association probabilities for the candidates, we 
compared radio positions and sizes available from Green’s catalog, with that of the GeV candidates’ 
localizations, localization errors, and extensions measured in this catalog. Using this set of spatial 
information, we derived two parameters. The first, Overlapjo^, provides a quantitative measure of 
whether the GeV localization is within the SNR’s angular extent in radio. This parameter, ranging 
from 0 to 1, is calculated as: 


Overlapio^. 


Radio n GeVioc 
mm{Radio, GeV\oc) 


(3.1) 
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Figure 1: Distribution of SNR radio diameters from Green’s catalog. The vertical dashed lines, the 
thin-dotted and thick-dashed lines indicate respectively the minimum detectable extension using 
lx, 10X and lOOx the isotropic background level, for more details see [4] and [1]. Roughly one 
third of the known Galactic SNRs may be resolved by the LAT if they are sufficiently bright GeV 
sources. 


where Radio represents the SNR’s radio disk and GeV\oc is the GeV 95% error circle. The no¬ 
tation X represents the area of X. The second parameter, Overlap^^t, quantifies whether the GeV 
candidate’s localization and extension is consistent with the location and extension of the radio 
SNR: 


Overlapext = 


Radio n G^Vext 
max{Radio^ G^Vext) 


(3.2) 


where Radio is again the SNR’s radio disk and G^Vext is the best fit GeV disk if the GeV candidate 
is significantly extended. If the GeV detection is consistent with a point source, we determined if 
the corresponding SNR’s radio size was consistent with it by redefining the Overlap^^t parameter 
as: 


Overlapext = 


Radio n G^Vnin 
Radio 


(3.3) 


where GeV^mn is the minimum resolvable radius, GeV^mn = 0.2° for this analysis, see also [4]. The 
values of Overlapi^^ and Overlap^^t for all significant GeV sources are shown in Fig. 2. We label the 
GeV detections with the most likely chance of true association as “classified candidates,” defined 
as those sources with Overlapext > 0.4 and Overlapj^e > 0.4. “Marginally classified candidates” 
are those GeV sources with a moderate chance of true association, defined as Overlapext >0.1 and 
Overlaptee >0.1 and at least one overlap estimator < 0.4. Candidates which do not have overlap 
parameters in these categories are referred to as “other” sources. 


4. Chance Coincidence 

In order to estimate the probability that any particular coincidence occurs by chance, we cre¬ 
ated a mock SNR catalog derived from Green’s catalog with SNR positions randomized excluding 
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Figure 2: Distribution of Overlapi^^ and Overlapg^t for all significant GeV sources. GeV candidates 
with significant extension are shown as open circles and GeV candidates consistent with the point 
hypothesis are filled. Stacked bar histograms of both parameters are also shown. The classified 
candidate region discussed is shown with a white background. Those points in the grey region to 
the upper right of the dotted line are marginally classified candidates. 

the original location and we run our standard pipeline and classification method. Any classified 
source found would be by chance. 

After the analysis, we find that only two out of 279 mock SNRs are spatially coincident with 
a GeV excess. Comparing the A^mock = 2 mock coincidences to the A/creen = 36 candidates passing 
the association probability we estimate a false discovery rate of ^ 6% for this specific realization 
of the mock catalog. We determined that, at 95% confidence, the number of false discoveries will 
be less than eight for any mock catalog prepared as described above, corresponding to an upper 
limit of 22% for the false discovery rate. 

5. Source of systematic errors 

For sources with significant emission, we estimated the systematic error propagating from the 
systematic uncertainty of the effective area and from the choice of lEM. For the former, we es¬ 
timated the systematic error by calculating bracketing IRFs, following the standard method [8]. 
While for the choice of the lEM we developped a new method. An earlier version of this method 
was described in [9]. Interstellar emission contributes substantially to EAT observations in the 
Galactic plane, where a majority of SNRs are located. Moreover, interstellar y-ray emission is 
highly structured on scales smaller than the regions of interest (Rols) typically used for this anal¬ 
ysis. To explore the systematic effects on SNRs’ fitted properties caused by interstellar emission 
modeling, we have developed a method employing eight alternative lEMs. The work in [10], using 
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the GALPROP^ cosmic ray (CR) propagation and interactions code, was the starting point for our 
alternative IBM building strategy. This strategy is different from what was adopted in the develop¬ 
ment of the standard IBM used in the rest of the analysis and in the usual Fermi BAT analysis. In 
order to build those eight models, we varied the values of three input parameters that were found 
to be the most relevant in modeling the Galactic plane: CR source distribution, height of the CR 
propagation halo, and Hi spin temperature [10]. The models were constructed to have separate 
templates for emission associated with gas traced by Hi and CO in four Galactocentric rings and 
an inverse Compton (IC) template covering the full sky. By allowing separate scaling factors for 
these different components of the model, we allowed many more degrees of freedom in fitting the 
diffuse emission to each Rol. 

Bor each significant candidate SNR we considered two hypotheses: the point source and one 
of the extended hypothesis that is preferred in the main pipeline analysis. Bor each significant 
candidate SNR we performed independent fits for each hypothesis for each of the eight alternative 
IBMs as well as the standard IBM, for a total of 18 fits of the each Rol. Bor each candidate SNR 
and alternative IBM we evaluate the preferred extension hypothesis, shown in figure 3a, and its 
significance, in figure 3b. 



Alternative lEMs with Significant Hypothesis Alternative lEMs with Extension Hypothesis Identical to Standard 

(a) Significance comparison (b) Best-fit model comparison 

Bigure 3: Bor each candidate’s analysis with each of the alternative IBMs, we computed the candi¬ 
date’s significance and the best extension hypothesis, either point-like or extended. Left: Number 
of analyses of each region with each of the alternative IBMs for which the candidate remained sig¬ 
nificant (Test Statistio 9). Right: Number of analyses of each region for which the best extension 
hypothesis with the alternative IBM remained the same as that found with the standard one. 


To identify which, if any, of these three sets of GALPROP parameters (CR source distribution, 
height of the CR propagation halo, and Hi spin temperature) has the largest impact on the fitted 
source parameters, for all the Rols we marginalized over the other two parameter sets and examined 
the ratio of the averaged candidates’ parameter values relative to the values’ dispersion. Bor a fitted 
source parameter P, such as flux, and a GALPROP input parameter set PaiEM = ^-g- spin 


^http://galprop.Stanford.edu/ 
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temperature Ts = {150 K, 10^ K}, the ratio is: 

\<Pi>-<Pj> 

R='—— - y- 

max[Gp.^apj) 


(5.1) 


In the previous definition < Pi > and Gp. are, respectively, the average and the standard deviation 
of the sample composed by the candidate parameters P obtained after the fit with GALPROP input 
models with the parameter i (e.g. spin temperature Ts = 150). A value of R >> 1 for a particolar 
GALPROP input parameter set means that all the other sets are negligible for that particolar can¬ 
didate. If in various Rols different GALPROP input parameter sets (CR source distribution, CR 
propagation halo height and Hi spin temperature) have a larger effect on P than the others, we need 
to test each source with all the alternative lEMs since in different part of the sky different maps are 
more relevant. As you can see in Figure 4a for the flux and 4b for the power-law (PL) index, none 
of the three GALPROP input parameters has R significantly smaller or larger than 1 such that the 
other GALPROP parameters can be neglected for all the candidates tested. 

We evaluated the systematic error due to the uncertainties of the interstellar medium modelling 
using the following formula. For each fitted parameter P we obtain a set of M = 8 values Pi that 
we compare to the value obtained with the standard model Pstd- Our estimate of the systematic 
uncertainty on P due to the modeling of interstellar emission is: 


/ 1 ^ 

t^sys,w — \ {Pj ~ PsTD )^• (5.2) 

V U i 

In Equation 5.2 the weights are w/ = 1 /af where (7/ is the statistical error of a parameter with a 
particular alternative lEM. 

We note that this strategy for estimating systematic uncertainty from interstellar emission mod¬ 
eling does not represent the complete range of systematics involved. In particular, we have tested 
only one alternative method for building the lEM and varied only three of the input parameters. 
This ensemble of models therefore cannot be expected to encompass the full uncertainty associated 
with the lEM. Further, as the alternative method differs from that used to create the standard lEM, 
the parameters estimated with the alternative lEM may not bracket the value determined using the 
standard lEM. Our estimate of the systematic error in Equation 5.2, accounts for this. 


6. Error comparisons 

The systematic error on flux estimated from the alternative lEMs tends to dominate over that 
due to effective area. This is true for most sources, regardless of their size, classification, or Galactic 
longitude. Moreover, the error propagated from the effective area tends to be < 10% of the flux 
for the large majority of candidates, while that from the alternative lEMs ranges over ^ 3 orders 
of magnitude, indicative of the diverse environments in which these candidates are found. We also 
observe that the candidates classified as other than SNRs, particularly when extended, tend to have 
somewhat larger alternative lEM errors than the classified and marginally classified candidates. 
For the PL spectral index the systematic errors from the bracketing IRFs dominate the alternative 
lEM errors for approximately half the candidates. Aside from candidates classified as other having 
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(a) Flux 
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T ▼ SNR G089.0+04.7 
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6 O SNR G347.3-00.5 
★★ Test SNRs’ Avg 
★★ All Candidates' Avg 
Marginal Avg 
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' Source distribution Halo height Spin temperature 


• • SNR G023.3-00.3 
▼ T SNR G089.0-H04.7 
< <1 SNRG120.1-t01.4 
► ► SNRG180.0-01.7 

• O SNR G213.3-00.4 
■ ■ SNRG260.4-03.4 
O O SNRG347.3-00.5 
irk Test SNRs' Avg 
irk All Candidates' Avg 
kk Marginal Avg 
irk Classified Avg 


(b)Index 


Figure 4: The impact on the fluxes and on the PL index of candidates for each of the alternative 
IBM input parameters (source distribution, halo height, and spin temperature), marginalized over 
the other GALPROP input parameters, is shown {R in the text). The stars represent the average 
ratio over the different candidates’ sets (classified, marginal, test and all). As no alternative IBM 
input parameter has a ratio significantly larger than 1, no input parameter dominates the systematic 
uncertainties of the fitted source parameter sufficiently to justify neglecting the others. 


somewhat larger systematic errors due to the alternative IBMs, systematic errors on neither flux nor 
index are significantly different among various subtypes: point or extended candidates, marginal 
or classified candidates, or young compared with interacting. 

The total systematic error, derived adding in quadrature the bracketing IRBs and alternative 
IBMs errors, dominates the statistical error on a candidate’s flux, as seen in Bigure 5a. The statisti¬ 
cal errors for PL index dominate over systematics for a number of candidates as shown in Bigure 5b. 


7. Conclusions 

We detected 102 candidates with a final source TS > 25 in the 279 SNR Rols (see Section 3). 
Of the 102 detected candidates, 36 passed the association probability threshold. Of these, 30 SNRs 
(^ 11% of the total) show significant emission for all alternative IBMs and are classified as likely 
GeV SNRs. An additional four were identified as sources which are not SNRs; two other can¬ 
didates were demoted to marginal due to their dependence on the IBM. Of the sources likely to 
be GeV SNRs, 17 show evidence for extension (TS^xr > 16). Only sources associated with SNRs 
G34.7—0.4 and GI89.1+3.0 show evidence of significant spectral curvature in the I — 100GeV 
range and are fit with logP spectra. Of the classified candidates, four extended and 10 point SNRs 
are new. Bor all the sources we evaluated their spectral and spatial characteristics with systematic 
and statistical uncertainties. A description of the catalog results in a multiwavelenght contest and 
their implication in the understanding of the CR acceleration can be found in [I] and [2]. 
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(a) Flux (b) PL Index. 


Figure 5: Comparison of the statistical and systematic errors, the latter derived from the alternative 
lEMs and the bracketing IRFs, for the flux and PL index. Open circles indicate extended SNRs 
while filled circles indicate point-like sources. All SNRs that passed classification are shown as 
black unless also classified as young non-thermal X-ray SNRs (blue) or as interacting with MCs 
(red). Candidates which did not pass classification but still had both fractional overlaps >0.1 are 
grey. If they are also young or interacting, they are outlined in blue or red, respectively, with the 
addition that all candidates classified as “other” are shown in green. For the flux the systematic 
error typically dominates the statistical error on the flux for all classes of candidates. For the index, 
in a number of cases the statistical error dominates the systematic. 
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